Quenching of magnetic excitations in single adsorbates at surfaces: Mn on 

CuN/Cu(100) 



Frederico D. Novaes 1 , Nicolas Lorente 2 , and Jean-Pierre Gauyacq 3,4 
1 Institut de Ciencia de Materials de Barcelona (CSIC), 
Campus de la UAB, E-08193 Bellaterra, Spain 
2 Centre d'Investigacio en Nanociencia i Nanotecnologia (CSIC-ICN), 
Campus de la UAB, E-08193 Bellaterra, Spain 
3 CNRS, Institut des Sciences Moleculaires d'Orsay, 
ISMO, Unite de Recherches CNRS-Universite Paris-Sud, 
Batiment 351, Universite Paris-Sud, 
91405 Orsay CEDEX, France 
4 Universite Paris-Sud, 
Institut des Sciences Moleculaires d'Orsay, ISMO, 
Unite de Recherches CNRS-Universite Paris-Sud, 
Batiment 351, Universite Paris-Sud, 
91405 Orsay CEDEX, France 

(Dated: June 30, 2010) 

The lifetimes of spin excitations of Mn adsorbates on CuN/Cu(100) are computed from first- 
principles. The theory is based on a strong-coupling T-matrix approach that evaluates the decay of 
a spin excitation due to electron-hole pair creation. Using a previously developed theory [Phys. Rev. 
Lett. 103, 176601 (2009) and Phys. Rev. B 81, 165423 (2010)], we compute the excitation rates 
by a tunneling current for all the Mn spin states. A rate equation approach permits us to simulate 
the experimental results by Loth and co-workers [Nat. Phys. 6, 340 (2010)] for large tunnelling 
currents, taking into account the finite population of excited states. Our simulations give us insight 
into the spin dynamics, in particular in the way polarized electrons can reveal the existence of an 
excited state population. In addition, it reveals that the excitation process occurs in a way very 
different from the deexcitation one. Indeed, while excitation by tunnelling electrons proceeds via 
the s and p electrons of the adsorbate, deexcitation mainly involves the d electrons. 

PACS numbers: 68.37.Ef, 72.10.-d, 73.23.-b, 72.25.-b 



I. INTRODUCTION 

Recently, a series of experimental studies^— with low- 
temperature STM (Scanning Tunnelling Microscope) re- 
vealed that isolated adsorbates on a surface could ex- 
hibit a magnetic structure, i.e. that a local spin could 
be attributed to the adsorbate. Interaction of this lo- 
cal spin with its environment results in several magnetic 
energy levels that correspond to different orientations 
of the local spin relative to the substrate. A magnetic 
field, B, was also applied to the system: increasing the B 
field decouples the local spin from its environment and 
the system switches to a Zeeman structure, thus help- 
ing to characterize the adsorbate magnetic anisotropy. 
The energy of the various magnetic levels were obtained 
via a low-T IETS (low- Temperature Inelastic Electron 
Tunnelling Spectroscopy) experiment, in which the tip- 
adsorbatc junction conductivity as a function of the STM 
bias exhibits steps at the magnetic excitation thresholds. 
The magnetic excitation energies are small, typically in 
the few meV range. Besides the spectroscopy properties, 
the IETS experiments also revealed that the conductance 
steps at the magnetic inelastic thresholds are very high, 
i.e. that the efficiency of the tunneling electrons in in- 
ducing magnetic transitions in the adsorbate is extremely 
large. In Mctal-Phthalocyanine adsorbates for exam- 



ple ■ , the tunneling current is dominated by its inelastic 
component when the bias is above the magnetic excita- 
tion thresholds. This efficiency is at variance with the 
efficiency of tunneling electrons in exciting vibrational 
modes in a molecular adsorbate which was observed and 
shown to only reach the few % ranged—. Several theo- 
retical accounts of the magnetic excitation process have 
been reported based on perturbation theory 2 *^— . or on 
a strong coupling approac h 13 ' 14 . In particular, the strong 
coupling approach quantitatively accounted for the very 
large efficiency of tunneling electrons in inducing mag- 
netic transition o 13 ' 14 . 

The existence on a surface of nano-magnets, the ori- 
entation of which could be changed at will by tunneling 
electrons, opens fascinating perspectives for the minia- 
turization of electronics. However, to lead to easily man- 
ageable devices, the excitation of local spins must have, 
among other properties, a sufficiently long lifetime. It is 
thus of paramount importance to know the decay rate 
of the excited levels of the local spin and in particular 
to decipher the various parameters and effects that gov- 
ern its magnitude. Experimentally, the local spins were 
observed in systems in which a coating on the surface 
was separating the magnetic adsorbate carrying the lo- 
cal spin from the metal substrate. Experiments on adsor- 
bates directly deposited on a metallic substrate did not 
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lead to sharp IETS structures^ as the others and this 
was attributed to a too short lifetime of the magnetic 
excitation on metals, stressing the importance of the de- 
coupling layer between local spin and substrate in sta- 
bilizing the magnetic excitation. Deexcitation of a local 
spin implies an energy transfer from the local spin to the 
substrate degrees of freedom, i.e. to the substrate elec- 
trons or to phonons. Phonons are not directly coupled to 
spin variables, but only via spin-orbit couplings (see e.g. 
a discussion in [Hj]). In contrast, the adsorbate spin vari- 
ables can be directly coupled to substrate electrons and 
electrons colliding on a magnetic adsorbate can easily in- 
duce magnetic transitions. Actually, this is exactly what 
happens in the magnetic excitation induced by tunneling 
electrons in the IETS experiments described above; in the 
de-excitation process the tunneling electrons are simply 
replaced by substrate electrons. The decay of excited 
magnetic states in individual adsorbates thus proceeds 
via electron-hole pair creation. Substrate electrons col- 
liding on the adsorbate can be thought to be as efficient 
in inducing magnetic transitions as tunneling electrons 
injected from an STM tip. In this qualitative view, one 
can expect the magnetic excitation decay rate to be the 
product of the collision rate of substrate electrons on the 
adsorbate by a very high efficiency factor. 

Recently, the decay rate of excited magnetic Mn atoms 
adsorbed on CuN/Cu(100) has been measured by Loth et 
alii via the analysis of the dependence of the adsorbate 
conductivity on the tunneling current. The decay of the 
magnetic excitations was interpreted in the above scheme 
as a decay induced by collision with substrate electrons. 
The lifetimes of magnetic excitations were typically found 
to be of the order of a fraction of ns. In the present paper, 
we report on a theoretical ab initio study of the lifetime 
of magnetic excitations in the Mn/CuN/Cu(100) system 
using both a DFT-based (Density Functional Theory) 
description of the system and the strong coupling for- 
malis m 13 ' 14 developed to treat magnetic transitions in- 
duced by tunneling electrons; the corresponding results 
are compared with Loth ct al data^i. 

II. METHOD 
A. Description of the magnetic deexcitation 

The present treatment of the decay of magnetic exci- 
tations closely parallels our earlier treatment of magnetic 
excitations in IETS (see details in Ref. [H[). We assume 
that the magnetic levels of the adsorbate can be described 
by the following magnetic anisotropy Hamiltonian^, 

H = g f i B B-S + DS 2 z +E(Sl-Sl), (1) 

where S is the local spin of the adsorbate, g the Landc 
factor and [Ib the Bohr magneton. B is an applied mag- 
netic field. D and E are two energy constants describing 
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FIG. 1: Energies of the six magnetic states in the 
Mn/CuN/Cu(100) system, eigenvalues of the magnetic 
anisotropy Hamiltonian (fl]). They are presented as a func- 
tion of the applied magnetic field, B, for a field along the 
x-direction. 



the interaction of S with the substrate, i.e. the magnetic 
anisotropy of the system. Diagonalisation of Hamilto- 
nian ([I) yields the various states, of the local spin 
of Mn (we have used S=2.5, g=1.98, D=-4I ,ueV and E=7 
/icV obtained in the experimental work-ii by adjustment 
to the magnetic excitation energy spectrum). For Mn on 
CuN, the local spin is 2.5 and there are thus 6 magnetic 
levels. 

Figure [T] presents the energies Ei of the anisotropy 
states, eigenstates of Hamiltonian (JTJ as a function of 
B, the applied magnetic field. Below, the ground state 
is noted and the excited states with i = 1-5. The easy 
axis of the system, the z-axis, is normal to the surface 
in this system. In the present study, coherently with the 
experimental study of Loth et al^, the B field has been 
put parallel to the surface, along the x-axis. As B is in- 
creased, the magnetic structure of the system changes 
from a magnetic anisotropy induced by the substrate 
with three doubly degenerate states at B=0 to the six 
states of a quasi-Zeeman structure at large B, where the 
states are eigenstates of S 2 and S x . The ground 
state at large finite B corresponds approximately to the 
M x = —2 state. The energy diagram in Fig. [T]thus cor- 
responds to the decoupling of the magnetic anisotropy 
by the B field. The structure appears a little complex at 
low B with several avoided crossings since the B field is 
not along the principal magnetic axis of the system. 

The aim of the present work is to compute the decay 
of the excited states by collision with substrate electrons, 
i.e. by electron-hole pair creation. The decay rate, Tx t,i, 
of an excited state, eigenstate of the Hamiltonian 

(P) with energy Ei, is the inverse of its lifetime n and 
it can then be written using matrix elements of the T 
transition matrix (we assume the energy variation of the 
T matrix to be small on the energy scale of the i — > f 
transition and we assume a vanishing temperature of the 
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substrate) 1 ^, 

1 ^ ^2^, 

/ / 

I _ 2 

x 2J \(kf,mf,&f\T\ki,mi,$i) 
k t ,kf 

mi ,m f 

x S(e l -e f )8(s l ~E F ), (2) 

where |<frj) are the final states of the decay, associated 
to an energy transfer of Sttf = E{ — Ef, and the total 
energy is Et = Ei + e% = Ef + ef. The initial and final 
states of the substrate electrons are noted by their wave 
numbers, fej and kf, and by their initial and final spin 
projections on the quantization axis, and m/. The 
substrate is assumed to be non-magnetic. Each term, 
Tij, in the sum over / is the partial decay rate of the 
initial state to a peculiar final state. 

The electron-adsorbate collisions being very fast, we 



can treat them without the magnetic anisotropy (Hamil- 
tonian ([!])) taken into account and later, include it in 
the sudden approximation. The electron-adsorbate col- 
lision is treated in a DFT-based approach (see section 
III Bl below) whereas the spin transitions are treated in 
the sudden approximation in a way similar to our ear- 
lier study of magnetic excitatio n 13 ! 14 . We thus define 
a collision amplitude for the substrate electrons inde- 
pendently of the magnetic anisotropy, it is a function 
of both the initial and final electron momenta, fcj and 
kf, and of the spin coupling between electron and ad- 
sorbate. The electron-adsorbate spin coupling scheme is 
defined via, St, the total spin of the electron-adsorbate 
system: St = S + s, where S is the adsorbatc spin and 
s is the tunneling electron spin. The projection of St 
on the quantization axis is noted My. If the adsorbate 
spin is S, there are two collision channels corresponding 
to St = S + \ and St = S — \ . The decay rate is then 
expressed as a function of the two channel amplitudes, 



1 v ■> 2ll8VLf X - 

f k } 



hi ,kf 
rrii ,77i f 



£ (kf,mf,$ f \ST,M T )T ST (S T ,M T \ki,m i ,$i) 

S t ,Mt 



6(Ei-e f )6(ei-E F ), (3) 



that can be re-expressed as 



r -E^F E S(e i -e f )6(s i -E F ) 



f 



ki,kf 
m. i ,m f 



^2{kf\T ST \k l )^2{mf,^f\ST,M T )(ST,M T \m l: <Pi) 

St A/It 



(4) 



One can see that the contributions from the two St 
terms are interfering. This corresponds to the a pri- 
ori general case where no selection rule apply to the 
tunnelling process. In that case, the above expresion 
can be simplified by making an extra statistical approx- 
imation that neglects the interferences between the two 
St channels. However, in our earlier studies on excita- 
tion processe s 13 ! 14 it was found that one of the two St 
channels was dominating the tunneling process between 
tip and adsorbatc. Here for deexcitation in the Mn on 



CuN/Cu(100) system, only one St coupling scheme does 
contribute significantly to the collision between substrate 
electrons and adsorbate (see Section III B|) and is thus in- 
cluded in the present treatment. One can note though 
that the dominating channel for the deexcitation process 
(collisions with electrons coming from the substrate) is 
different from that dominating the excitation (collision 
with tunnelling electrons) (see Section III Bp . 

If one considers a single St channel, the decay rate can 
be rewritten as: 



J 



1 _ ^ 27T<ft2/ 



E 

/ 



h 

2n5VL f 



^2{mf,^f\S T ,M T ){S T ,M T \mi,^i 



Mi 



f k i: k f 
mi ,m / 

S( ei -e f )5( ei -E F ) \(kf\T s -\h)\ 2 Pf^ Lmf , 



ki,kf 
rrii,m f 



(5) 
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where Pfl. 



>/ m * s a ma t rrx element of spin variables only, corresponding to the St coupling scheme for the (i, rrii 



f,mf) transition. Equation (|5J then becomes: 



f 



= E^ I *w| E ^ 



The total and partial decay rates are then obtained via: 




Ttii — >/,m f 



(6) 



^ = E r W = E ^ £ (2 7 r) 2 T s -(^)P Spi „(5 T ,i 



(7) 



Each partial decay rate thus appears as the product of 
a spin transition probability by the electron flux hitting 
the adsorbate in the energy range able to perform the 
studied transition. 

As we will see below, the DFT-based calculation of 
the electron flux is not performed in the St, Mt spin 
coupling base, but by specifying the scattering electron 
spin state (majority or minority). Thus we cannot per- 
form a one to one identification. However, scattering 
through the adsorbate being dominated by a St chan- 
nel , we can identify the electron flux in the St channel 
(T St (E f )) with its equivalent in the DFT-approach, the 
total electron flux hitting the adsorbate T Total (E F ). So 
that finally, the total decay rate is obtained as: 



1 

Ti 



= Er W = T™«\E F ) E S -^Ps P in(i -> /)(8) 



The decay rate is then equal to the total flux of sub- 
strate electrons hitting the adsorbate per second in the 
appropriate energy range, T Total (E F )5£lf, times a spin 
transition probability. Below, T Total (EF) is identified 
with the equivalent quantity computed in the DFT ap- 
proach (Section lllBp . One can stress the great similarity 
of this expression with that of the inelastic conductivity 
obtained in Ref. [l3lfl4j as the total conductivity times 
a spin transition probability. 



DFT-based calculation of the substrate electron 
collision rate 



T{E F ) = (2tt) 



E 

ki,kf 



(ki\f\k f ) 6(ei-ef)6(ei-E F ) (9) 



In this equation, \ki) and \kf) arc asymptotic states that 
are solutions to the Hamiltonian sufficiently far away 
from a scattering center. In our case, the scattering 
center is the adsorbed atom. The goal is to calculate 
the amount of substrate (CuN/Cu(100) surface) electrons 
that scatter at the adsorbed Mn atom. From this quan- 
tity, we obtain T Total (E F ) by summing the values for 
both spins, 



T Total {E F ) = T(E F ,t) + T(E F ,l) 



(10) 



Following standard scattering theory formalism, we 
write the solutions to the full Hamiltonian (sur- 
face+adsorbed atom) as, 



|V) = \ki) + G r V\h) = \h) + G 0r V\tP) (11) 

where V is the perturbation that couples the adsorbed 
atom and the substrate; and G r and G 0r are the full sys- 
tem's and unperturbed retarded Green's functions, re- 
spectively. The T operator is then written as, 



T = V + VG r V 



(12) 



In order to actually compute these vectors, matrices, 
and ultimately T(E F ), we have used the SIESTA^ code, 
and a set of strictly localized atomic-like orbitals as a 
basis set. In terms of these, the Hamiltonian has the 
form, 



The quantity that we want to compute from first- 
principles simulations is, 



H 



Has H^a 



(13) 
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where the H M/J are themselves matrices whose elements 



H 



DFT 



jj — are computed with orbitals cen- 

tered around either surface (S) atoms or around the 
adsorbed atom (A). These orbitals are not orthogonal, 
(4>i\4>j) = Si j. In this representation, the coefficients of 
the and \ki) states form a vector of components^, 



CA 



u 5 





(14) 



(15) 



G 0r (E) = (E + S c 



Hn 



E + = lim E 
8->o 



(21) 



(22) 



What we need next, is to determine the perturba- 
tion matrix V entering equation (|11| . By multiplying 
equation (fT9f by (E + S — H) from the left, we get that 
Vss = and Vas = H^s — ESas, where we have used 
equation (fT6|) and the fact that, 



and the are solutions to, 



lim idScfr. 







(23) 



(Ho-£So)c fe! =0 



Hn 



So 



H ss 




Sss 
S AA 



The matrix version of equation is, 



CA 



^SS ^SA 

Gr c*r 

AS ^AA 

Vss Vsa 

V AS Vaa 



(16) 



(17) 



(18) 



c s 




(19) 



and the Green's functions matrices are defined as, 



Using the second part of equation pip - the one that 
involves G 0r - in a matrix form equivalent to equation 
(fl"9|) . we get that Vaa = 0. Imposing V to be hermitian, 
we finally have, 



V = 



(H 



AS 



ES 



AS) 



(Hsa — ESsa) 




(24) 



and we can see that (Ho + V — ESo)c^ = 0. The energy- 
dependent form of V can be traced back to the the use 
of a non-orthogonal basis se t 22 ' 23 . 

The value of (ki\T\kf) can thus be calculated with, 



(ki\T\k f ) 



( T ss 


T S a \ 


(1) 


{ T AS 


Taa J 





c s T ssc s 

4**(V S AG r AA V A s)c k s f 



(25) 



G r (E) = (E+S-H) 1 



where we have used equation (|12[) . Substituting this back 
(20) into equation ([9]), we get that, 



= Y, c s*( Y sAG A AVAs)4ulr(VsAG AA VAs)4>5(e l -e f )6(s l -E F ) 1 

ki,kf 



= Tr 



Vas (Ylc^c^Siei-E^ VsaG aa Vas c*'c*'**( £< - e f )j V SA G AA 



(27T) 



-Tr [T AA {E F )G AA {E F )T AA {e l = E F )G a AA {E F ) 



and we obtain the final expression for T(E F ) 



T{E F ) = Tr [r AA (E F )G r AA (E F )r AA (E F )G a AA (E F ) 



(26) 



(27) 



In the derivation of equation (f2l))) we have used that: 



1. CgCg* can be considered to define a matrix: Cg be- 



6 



ing a column vector, and Cg* a row vector. 

2. The cyclic properties when taking the trace of a 
product of matrices. 

3. The discontinuity of the retarded and advanced 
Green's functions at the real axis 2 ^ that in our basis 
set, considering item 1, gives, 

44'8(E - e 3 ) = ^(G%(E) - G s a s (E)) (28) 

3 

4. The self-energy and gamma matrices are written 
as2£, 

V AS G° s f ] V SA = T?J$ (29) 

T AA = i [E r AA - £^ A ] (30) 

We now describe the procedure to obtain T(E F ), de- 
fined in equation (|27p . from our ab-initio simulations. 
The ideas are the same as the ones used to derive the 
equations for the transmission function for electronic 
transport calculations, but since the final formulas are 
not exactly the same, we here include the derivation of 
the relations that we have used in the present work. To 
do so, we start by re- writing the DFT Hamiltonian (equa- 
tion dm), 

(~H.ee Hec \ 
H CE H cc U CA (31) 
H A c Haa J 

where Hee is a semi- infinite matrix describing the semi- 
infinite electrode (a region where the electronic structure, 
and matrix elements are assumed to be already bulk- 
like); and Hcc is the "slab" that describes the actual 
surface. The size of the C (for contact) region is in prin- 
ciple arbitrary (but finite!), as long as it is thick enough 
to have: i) the He A (H AF ) matrix elements equal to 
zero; and ii) the Hee matrix elements sufficiently con- 
verged to bulk values. 

For the purpose of simplifying the notation, let us de- 
fine an h matrix to be, 

h = ES - H (32) 

To compute T(E F ), we need G 7 AA (E F ) and T AA {E F ), 
where, 

G r AA = (h AA - Haa) (33) 

From Eqs. (|2"9"|) . (|3"0"|) and (|B"3"|) . we can see that what we 

need is to compute G'qq the finite portion of G°aa 
needed to compute the self energies, 



G°cc = (hcc — hcsG^gjhec) (34) 

where it is important to note that only a finite number of 
elements of hcE (h_ec) have non-zero values, hence only 
a finite portion of G E r E needs to be calculated. Since the 
\\ee matrix elements are assumed to be already "bulk- 
like" , the G E r E matrix is obtained using the matrices 
extracted from a bulk calculation of the electrod o 25 ' 26 . 

From these equations, we see that T(E F ) represents 
the flux of electrons coming from the substrate and scat- 
tering off the adsorbate back into the substrate. This 
quantity is different from the transmission function ap- 
pearing in a Landauer-like approach 2 ^, where the adsor- 
bate is connected to two electrodes. The difference is 
clear in the above equations, here a unique reservoir (the 
substrate) appears in the self-energies, and hence the de- 
cay appearing in the elastic Green's function, Eq. (|33p . 
has only one self-energy instead of two in the transport 
case 2 ^. 



III. DENSITY-FUNCTIONAL STUDY OF A 
SINGLE Mn ADSORBATE ON A CuN/Cu(100) 
SURFACE 

The ground-state electronic-structure configuration 
and the value for T Total (E F ) were obtained by density- 
functional-theory (DFT) simulations. Our DFT calcu- 
lations were performed using the SIESTA code2£. The 
super cell contained at least six 4x4 Cu (100) layers - 
a larger number of layers was used to test the conver- 
gence with respect to the size of the C region, as dis- 
cussed in section Hi Bl one CuN layer; and one Mn atom 
(see Fig. The Mn atom and the two outermost lay- 
ers were relaxed until the forces were smaller than 0.03 
eV/Ang. A sampling of 3 x 3 x 1 k-points was used. We 
have used the generalized gradient approximation 2 ! for 
the exchange-correlation potential. 

We have obtained a total spin polarization (Q up — 
Qdown) of 4.7 Bohr magnetons, which essentially corre- 
sponds to a S = | spin configuration, in good agreement 
with experiment^. The spin polarization is localized 
around the Mn atom, that has five half-filled d-orbitals, 
corresponding to five unpaired electrons. This can clearly 
be seen in Fig. [3] 

We now turn to the calculation of T Total (E F ). As 
already discussed in section IIIB| the thickness of the C 
region of Eq. pip is arbitrary (to some extent), and here 
we have taken the first three surface layers (including the 
CuN layer), as shown in Fig. [2] Since there is also some 
arbitrariness in the definition of what the atom is, we 
have considered different basis sets: from SZP to DZP, 
with larger and shorter cutoff radii. The differences in 
the results were not appreciable, and so here we report 
just the results for a DZP basis with large cutoff values 
— up to 10 A for the Mn basis orbitals. 
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FIG. 2: Top and lateral view of the geometry used to describe 
the Mn adsorbate (green) on top of a CuN/Cu(100) surface. 
The Mn atom is located on top of a Cu atom of the CuN 
layer, pushing it down while pulling up the N (blue) atoms 
close to it. 
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FIG. 3: Projected density of states (PDOS) on the Mn atomic 
orbital for the case of Mn on CuN/Cu(100). For all the curves 
shown here, the positive (black) curves corresponds to the 
majority spin, and the negative (red) to the minority spin. 



We have computed T(E) for a range of energy values 
around the Fermi level, and the results are show in Fig. |4] 
As for T Total (E F ), we have found T Total (E F ) = T(Ep, t 
) + T(E F , i) = 0.8 + 0.3 = 1.1. When we analyze the or- 
bital contribution to the function T(E) of Fig. [4] we iden- 
tify the d xz orbital of Mn as responsible for the electron 
scattering off the Mn adsorbate. The above information 
on the characteristics of the substrate electrons hitting 
the adsorbate can be used to further specify the inputs 
of our spin-transition calculations. Indeed, the ground 
state of the system corresponds to putting one electron 
into each Mn d orbital (defined with the appropriate sym- 
metry). This generates the S = 5/2 M = 5/2 state of 
the adsorbate. The dominant contribution to substrate 
electrons going through the adsorbate is found to involve 
the d xz orbital with majority spin as the transition in- 
termediate; this is interpreted as a process involving a 
positive ion intermediate of St = 2 symmetry. Simi- 
larly, the contribution associated to substrate electrons 
going through the d xz orbital with minority spin is inter- 
preted as a process involving a negative ion intermediate 
of St = 2 symmetry. So in all cases, the deexcitation 
process induced by substrate electrons going through the 
Mn adsorbate involves a St = 2 intermediate and the 
associated electron flux is given by T Total (Ep) = 1.1. 

At this point, one can stress the stark contrast between 
the present deexcitation study and our earlier study on 
magnetic excitation by electron tunnelling between the 
tip and the substrate (Rcf. fl3l]). In the tunnelling elec- 
tron case, the Mn orbitals contributing to the transmis- 
sion were the extended s and p orbitals whereas here, for 
the electrons scattering from the substrate into the sub- 
strate via the adsorbate, a d orbital is dominating. In 
addition, the spin symmetry of the scattering intermedi- 
ates are different : St — 2 vs St = 3. This study permits 
us to conclude that the de-excitation of spin states via 
electron-hole pairs takes place through the Mn d elec- 
trons and in particular the d xz orbital, while the spin 
excitation process proceeds via the tunneling electrons 
and are of s and p characters. 



IV. LIFETIME OF THE EXCITED MAGNETIC 
STATES 

Equation ||HJ) together with the results of section lllll has 
been used to compute the decay rates of the excited states 
of the system. Figure [5] presents the total decay rate, 
TTot,i, inverse of the lifetime, of the five excited states as 
a function of the applied magnetic field B. Two different 
regimes can be seen: low B and large B. At large B, the 
magnetic structure of the system is quasi-Zeeman. As a 
consequence, the decay of an excited state is dominated 
by one channel corresponding to a AM X = — 1 selection 
rule. The B dependence of the decay rate is then that 
of the energy change associated to the decay, i.e. it is 
linear in B with a slope proportional to the dominant spin 
transition probability, Ps P m- At small B, the variation 
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FIG. 4: T(E) for the case of Mn on CuN/Cu(100). The 
black (dashed) curve corresponds to the majority spin, and 
the red (continuous) to the minority spin. T Total (E F ) is the 
sum of T(E) for both spins at E = Ef- A sampling of 15 x 15 
k-points was used to obtain converged values for T(E). 
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FIG. 5: Decay rate (in meV) of the five excited magnetic 
states in the Mn/CuN/Cu(100) system as a function of the 
applied magnetic field B. The B field is along the x-axis. The 
decay rate is the inverse of the excited state lifetime. 



is more complex, reflecting the complex decoupling of 
the anisotropy by the B field on the x-axis (see Fig. [TJ . 
However, one can notice that the lowest excited state 
remains quasi-degenerate with the ground state almost 
up to IT (Fig.[lJ, so that its decay rate is extremely small 
due to a quasi- vanishing 60. f (see Eq. ©). Actually, this 
simply means that at low B, for a finite temperature, the 
two lowest states arc roughly equally populated. As for 
the states 2-5, at low B, their decay rate is in the 2.0 [ieV 
range, corresponding to a lifetime of the order of 0.3 ns. 

One can stress that a different direction of the B-field 
leads to a different behaviour of the decay rates. The 
limit B = is the same, obviously and the large B is 
almost independent of the B direction in the present case, 
but not completely due to an incomplete decoupling of 




FIG. 6: Comparison between the present decay rate 
(full lines) of the five excited magnetic states in the 
Mn/CuN/Cu(100) system with the experimental results of 
Loth et a}i£ (symbols). The black symbols and line corre- 
spond to an applied B-field of 7T and the red symbols and 
line to 3T (the B-field is along the x-axis). For the sake of 
comparison, the experimental results of Loth et al have been 
multiplied by a global factor equal to 3.1. 



the magnetic anisotropy for B = 7 T. The decay rates 
in the intermediate region, where the decoupling of the 
anisotropy occurs, depend of the direction of the B field 
with respect to the magnetic axis of the system. 

Figure [6] presents a direct comparison between the 
present decay rate of the five excited states of Mn on 
CuN/Cu(100) and the decay rate extracted by Loth et 
aU£ from their experimental data. Two values of the 
magnetic field are presented: 3 and 7 T. For the sake of 
comparison, the experimental results of Loth et al have 
been multiplied by a global factor equal to 3.1. It ap- 
pears that the present study reproduces extremely well 
the state dependence and the magnetic field dependence 
of the excited state lifetimes. However, the present re- 
sults for the decay rates are a factor 3 larger than the 
experimentally extracted data. Besides inaccuracies and 
approximations in the experimental and theoretical pro- 
cedures, one can invoke the sensitivity of the present re- 
sults on the energy position of the d xz orbital in the cal- 
culations. Indeed, as we can see in Fig. 3 and 4, the 
electron flux hitting the adsorbate at Fermi level corre- 
sponds to the tail of the d xz majority and minority spin 
orbitals and any inaccuracy in the orbital energy directly 
affects, T(Ep), the electron flux. 



V. MODELLING OF THE Mn JUNCTION 
CONDUCTANCE AT LARGE CURRENT 

The experiments of Loth et alii introduce two main 
effects compared to the earlier ones in Ref.Q. First, by 
using a polarized tip, these experiments introduce an un- 
balance between the spin up and spin down tunneling 
electrons that reveals the spin-dependence of the junc- 
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tion conductance. Second, they consider large currents 
flowing through the junction allowing a stationary pop- 
ulation of excited states induced by the tunneling elec- 
trons. Indeed the two effects are linked together, tunnel- 
ing electrons of different spins leading to different excited 
state populations. Below, we examine these two effects 
in the Mn on CuN/Cu(100) system allowing a modelling 
of the experimental situation. 



A. Spin-dependence of the conductance of the 
adsorbed Mn atom 

In references [lolflij . we developed a strong coupling 
treatment of the elastic and inelastic tunneling of elec- 
trons through a magnetic adsorbate that inspired the 
above treatment of excited magnetic state decay. As ex- 
plained in section UH the magnetic anisotropy terms are 
small, and one can treat them in the sudden approxi- 
mation and define a tunneling amplitude, independent of 
the magnetic anisotropy. As the main result (see details 
in [lj]), the system conductance for the Mn atom in the 
magnetic state i, (in the case of only one St symmetry 



1.5 



contributing to tunneling) is given by: 



dJ 
dV 



V A ■ ■ A* 

Z^ij ■ n J,h">' Tl j,n,m' 



(35) 

where Co is a global conductance; on the small energy 
range that we consider here and for a fixed tip-adsorbate 
distance, Co can be considered as constant. V is the 
junction bias and E n the energy of the various magnetic 
states, n, of the system. The Aj >n _ m coefficients are spin 
coupling coefficients giving the expansion of the eigen- 
states of S T and St,z on the initial or final states of 
the tunneling process, \m,<f> n ), (m refers to the tunnel- 
ing electron spin state and tfi n is the magnetic anisotropy 
state of the adsorbate, eigenstate of the Hamiltonian (JTJ): 



A, 



{S T ,M T \m,<j) n ), with j = (St, Mi 



(36) 



In both treatments (inelastic tunneling and excited state 
decay), the final result appears as a product of a 
magnetism-free quantity by a spin-coupling coefficient 
term, i.e. a global conductivity of the system is shared 
among the various anisotropy channels. The excitation 
probability is then only dependent on the weight of the 
incident and final channels in the intermediate tunneling 
intermediate and it can be very large. Such an efficient 
inelastic process has been invoked in several other pro- 
cesses involving angular momentum transfer (rotational 
or spin) in gas phase or surface problems^— , in all 
cases, it lead to high probabilities of inelastic scatter- 
ing (sec also a discussion in Ref. [HI). A change in 
the adsorbate-tip distance leads to a change in Co only 
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FIG. 7: Conductance of the various magnetic states in the 
Mn/CuN/Cu(100) system as a function of the tip bias in mV 
for a fully polarized electrode (spin down). The B-field is 
along the x-axis and is equal to 3T. The conductance has 
been normalized in such a way that the conductance for the 
ground state and a non-polarized tip is equal to 1 at zero bias. 



and consequently in the tunneling current, but without 
a change in the excitation probabilities. 

In the case of Mn on CuN/Cu(100), a DFT calcula- 
tion showed that the St = 3 symmetry is dominating 
the tunneling process^ and this accounted well for the 
observations of Hirjibchcdin et al^, obtained with non- 
polarised tunneling electrons. In Ref. (l3lfl4j . we only 
considered tunneling of non-polarised electrons, i.e. we 
summed the contributions from the two electron spin di- 
rections, both in the incident and final channels (the sum 
over m and m' in equation (|35|l ). Here we consider tun- 
neling for a fixed direction of the electron spin (fixed 
m or m') in the incident or final state (the spin direc- 
tions are defined along the x-axis parallel to the applied 
B field). We thus use equation (|35|) with the sum over 
77i (or the sum over m') removed from the numerator. 
Figures [7] and [8] present the relative conductance of the 
various magnetic states of Mn for a fully polarized elec- 
trode and a B field of 3 T. Figure [5] corresponds for V > 
to an incident electron with an 'up' spin and for V < 
to an electron in a final 'up' state. Fig. [7] presents the 
'down' equivalent. In both cases the conductance is nor- 
malised in such a way that the conductance at V = 
for a non-polarised beam is equal to 1. The conductance 
presents steps at the magnetic excitation thresholds, due 
to the excitation induced by the tunneling electrons; the 
conductance also takes into account the possibility of de- 
excitation processes induced by the tunneling electrons, 
these present no energy thresholds. No broadening effect 
has been introduced in the conductance in Figs. I7|8I9I 
and the inelastic steps should be vertical; the finite slope 
visible in the figures comes from the finite number of V 
points that were actually computed. 

It appears that the conductances of the various 
anisotropy states for the two spin directions are quite dif- 
ferent both in magnitude and shape (some have inelastic 
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FIG. 8: Conductance of the various magnetic states in the 
Mn/CuN/Cu(100) system as a function of the tip bias in mV 
for a fully polarized electrode (spin up). The B-field is along 
the x-axis and is equal to 3T. The conductance has been nor- 
malized in such a way that the conductance for the ground 
state and a non-polarized tip is equal to 1 at zero bias. 
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FIG. 9: Conductance of the various magnetic states in the 
Mn/CuN/Cu(100) system as a function of the tip bias in mV 
for a non-polarized electrode. The B-field is along the x-axis 
and is equal to 3T. The conductance has been normalized in 
such a way that the conductance for the ground state and a 
non-polarized tip is equal to 1 at zero bias. 



contributions and some do not). This is not surprising 
within our strong coupling approach: the magnitude of 
the conductivity is directly given by the weights of the 
initial and final channels (the A^ n ^ m coefficients) in the 
tunneling symmetry [St — 3 in the present case) and 
these are strongly dependent on the considered initial 
and final states. One can also notice that the conductiv- 
ity shown in figures [7J and [5] are discontinuous at V = 
(except for the ground state) due to the switch of def- 
inition at V = (spin selection for the initial state vs 
spin selection for the final state). The behaviour seen in 
[7J and H] can be easily understood for the 'large' B case 
depicted there. In that case, the anisotropy states are 
roughly Zeeman states, eigenstates of S x with the eigen- 



values M x . In that case, the increasing order of energy of 
the magnetic states corresponds to the increasing order 
of M x and the transitions induced by a tunneling elec- 
tron verify the AM X = ±1 selection rule (the rule is strict 
only in the perfect Zeeman limit). As a consequence, ex- 
citation corresponds to a AM X = +1 selection rule and 
it can only exist for an incident 'up' electron and an out- 
going down electron, and this appears clearly in Fig. [7J 
and[H where the conductivity exhibits an inelastic step 
for only one sign of V, different for spins up and down. 
The inelastic steps appear as different energies for the 
different excited states, this is due to the fact that at 3 T 
the structure is not yet a perfect Zeeman structure (see 
Fig. [1]), in this limit all inelastic steps would be at the 
same position given by gfisB. Similarly, the very small 
excitation steps appearing at higher energy in addition 
to the AM X = +1 selection rule steps, as well as those 
appearing in the 'forbidden' V side, are due to the small 
difference from a pure Zeeman structure. The disconti- 
nuity in the conductance at V = is due to the existence 
of de-excitation processes induced by the tunneling elec- 
trons. Similarly to the excitation processes (well visible 
in Fig. [7Jand[5]), these arc highly dependent on the sign of 
V, leading to a discontinuity at V = 0. The purely elas- 
tic conductance is continuous at V = 0, as is the ground 
state conductance. Figure [TU] shows a qualitative picture 
explaining this bias asymmetry. 

Figure [S] presents the relative conductance for the var- 
ious magnetic anisotropy states in the case of a non- 
polarised tip. Similarly to Fig. [Jj and [H the conduc- 
tance has been normalised so that the conductance at 
V = for the ground state is equal to 1. All the con- 
ductivities are now continuous at V = and symmetric 
in V since no selection rule is imposed on the tunneling 
electron spin. The conductivity for each magnetic state 
exhibits a large inelastic step (except the highest lying 
state!) corresponding to the AM, = ±1 selection rule; 
much smaller excitation steps also appear at higher ener- 
gies (barely visible in the figure) due to the non-perfect 
Zeeman structure. For all excited states, inelastic tun- 
neling is significant compared to elastic tunneling, it is in 
the 15-25 % range, and different for the different states. 
One can also notice in Fig. |2]that the conductivities at 
large bias for all the states are equal, as a consequence 
of the summation over all possible excited states in Eq. 
1331 One can then conclude from Figs. OH] that the con- 
ductivity in the case of a polarised tunneling electron 
is quite different from that for a non-polarised one; in 
general, simply looking at majority and minority spins is 
insufficient, one has to resort to spin coupling arguments. 



B. Description of the conductivity in the presence 
of excited states 



The excitation of magnetic states by tunneling elec- 
trons is very efficient in the Mn/CuN/Cu(100). Though 
not as efficient as in some metal-phthalocynanine 
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FIG. 10: Qualitative scheme of the bias asymmetry in the tunneling process. For clarity of the illustration, we assume a, S — ^ 
adatom in its 'spin-up' state. The tip is completely polarized along the 'up' direction, (a) For positive bias, V > 0, electron 
tunneling can only be elastic because the tip is spin polarized and aligned with the adatom's spin, (b) For negative bias, V < 0, 
unpolarized electrons flow from the substrate and they can tunnel either elastically or inelastically. Hence, tunneling is different 
in cases (a) and (b) and so it presents a bias asymmetry due to the spin polarization of the tip. 



case^^, this can lead to a significant stationary pop- 
ulation of excited states in an experiment with a finite 
current. This effect has been very clearly demonstrated 
by Loth et alii. We modelled this situation in a way 
very similar to that used by Loth et ali£. The main dif- 
ference lies in the absence of adjustment parameters in 
the present work: with the present treatment of the ex- 
citation process (section IV A[) and of the decay process 
(sections HI1 and HV|) . we can quantitatively predict the 
behaviour of the conductance as a function of the tun- 
neling current. 

The basic idea is to compute the stationary popula- 
tion of excited states that is induced by the tunneling 
current and then, via the excited conductance discussed 
in IV Al to get the junction conductivity. For an STM tip 



positioned above the Mn adsorbate, a tip bias V and a 
tunneling current I, the time dependence of the popula- 
tion, Pj (J, V), of the magnetic state, i, is given by: 



+ ^/';Mi ;/;,(/. li • (37) 

3 

where Fjj is the partial decay rate of state i towards state 
j (sections HTl and HV|) . F i: j(I,V) is the transition rate 
from state i to state j induced by the tunneling electrons 
(section [VA| . It is given by: 
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where Pspinihj) is the spin coupling coefficient of the 
considered i — » j transition (section IV A[) for the con- 
sidered spin states of the tunneling electron. Aij(V) is 
an energetic factor; for the vanishing temperature consid- 
ered here, it is equal to V + Ei — Ej for an open excitation 
channel (i — > j transition), to for a closed excitation 
channel and to V for a de-excitation channel. C is a factor 
corresponding to the global conductance of the system. 
Below, the factor C is set so that the junction conduc- 
tance at V = 0, G, has a fixed value and then the whole 
spectrum of conductivity as a function of V is computed; 
changing C corresponds to moving the tip with respect 
to the Mn adsorbate i.e. to changing the magnitude of 
the tunneling current. 

Equation (|37|) yields the time dependent population 
of the excited states, including the transient regime at 
the switching of the applied bias. The experiments be- 
ing slow on the time scale of the excited state relaxation 
time (typically a fraction of ns, as seen in Fig. [5]), the 
populations quickly reach stationary values which deter- 
mine the observed conductivity. These stationary values 
are obtained by solving the homogeneous set of equa- 
tions obtained from (|37|) by setting all time derivatives 
to zero. Once the populations are known the junction 
conductance is obtained by summing the contributions 
of the different states. 



C. Population of the excited states 

Loth et al^ performed their experiments with a partial 
polarisation of the tip typically equal to 77 = 0.24 so that 
the two spin directions of the electron have probabilities 
equal to 0.5 (1 ± 77). In the present system, the ground 
state at large B is almost the M x = — § state. The 
polarisation of the tip is in the same direction, so that 
electrons with spin down are dominating at V > and 
holes with spin down at V < 0. 

The corresponding population of excited states for a 
junction conductance at V — equal to 2.10 -6 S, a 
magnetic B field of 3 T and a tip polarisation equal to 
rj = 0.24 are shown in Fig. [TT] as a function of V. Be- 
low the threshold for the — )• 1 transition, the popula- 
tion is entirely in the ground state, beyond this threshold 
the excited state population quickly increases as |V| in- 
creases. One can notice that the excitation process at 
3T basically goes step by step from state to 5 follow- 
ing the energy (and index) order, so that excitation of 
the higher lying states is a multiple order process and 
rises much more slowly than excitation to state 1. At 
large |V|, the system reaches a regime where the popula- 
tion is independent of V, it corresponds to the situation 
where transitions (excitation and de-excitation) induced 
by tunneling electrons dominate over the excited state 
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FIG. 11: Population of the six magnetic states in the 
Mn/CuN/Cu(100) system as a function of the tip bias. The 
magnetic field is along the x-axis and equal to 3T. The tip po- 
larisation, rj, is equal to 0.24. The conductance at zero bias 
is equal to 2.10" 6 S. 



decay. Excitation by electrons and holes are very differ- 
ent leading to very different excited state populations in 
the V > and V < cases. 'Down' polarised holes have 
a much higher excitation efficiency than 'down' polarised 
electrons, and this is a direct consequence of the effects 
observed in Fig. [JJ and [5] and discussed above. Actually, 
if we forget the very tiny excitations that correspond to 
the non-perfect Zeeman limit, a fully polarised tip would 
only lead to an excited state population for V < (see 
Fig. [7]); here with a finite 77 polarisation, the excited 
state populations at V > are due to the minority spin 
direction electrons in the tip. Indeed for a non-polarized 
tip, excitation by electrons and holes are equivalent, con- 
sistently with Fig. [SJ 



D. Modelling of Loth et al experiment 

Figures [JJ] and [T3] present the conductivity of the Mn 
junction computed with a finite excited state popula- 
tion for different values of the conductance at V = : 
0.1, 0.2, 0.5,1., 2., 5. and 10. (all in 10" 6 S). They corre- 
spond to a B field of 3 and 7 T, resp. and a tip polari- 
sation of 0.24. The computed conductivity as a function 
of V has been convoluted with a Gaussian of 0.2 meV 
width to mimic the various broadening phenomena. Note 
that the actual height of the inelastic conductivity step 
for B = 3 T depends on the magnitude of the broaden- 
ing effect, due to the small excitation energies (Fig. [TJ). 
The behaviour is the same as observed experimentally: i) 
there is a sharp peak at very low voltage corresponding 
to the excitation thresholds; ii) the conductivity drops 
as |V| is increased due to the excited state population; 
Hi) this drop is steeper and deeper on the V < side; 
iv) the steepness of the conductivity drop decreases as 
the tunneling current decreases to practically vanish for 
very small conductances and v) the conductance drop is 
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FIG. 12: Relative conductance of the Mn/CuN/Cu(100) sys- 
tem as a function the tip bias. The tip polarization is r\ — 0.24 
and the B field, equal to 3T, is along the x-axis. The finite 
population of the excited states is taken into account. The 
various curves correspond to various absolute conductances 
at zero bias (0.1, 0.2, 0.5, 1., 2., 5. and 10. 10" 6 S). In the 
figure, the conductance is plotted in relative value, with the 
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FIG. 14: Relative conductance of the Mn/CuN/Cu(100) sys- 
tem as a function the tip bias. The tip is non-polarized and 
the B field, equal to 3T, is along the x-axis. The finite pop- 
ulation of the excited states is taken into account. The con- 
ductance is independent of the current. In the figure, the con- 
ductance is plotted in relative value, with the conductance for 
zero bias set to 1. 
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FIG. 13: Relative conductance of the Mn/CuN/Cu(100) sys- 
tem as a function the tip bias. The tip polarization is r\ — 0.24 
and the B field, equal to 7T, is along the x-axis. The finite 
population of the excited states is taken into account. The 
various curves correspond to various absolute conductances 
at zero bias (0.1, 0.2, 0.5, 1., 2., 5. and 10. 10" 6 S). In the 
figure, the conductance is plotted in relative value, with the 
conductance for zero bias set to 1. 



steeper at 3T than at 7 T due to a longer lifetime of 
the excited states at 3T. The height of the inelastic steps 
and their asymmetry around zero, as well as the extent 
of the decrease of the conductance as \V\ increases de- 
pends on the polarisation of the tip. In the case of an 
unpolarised tip (?y = 0), Fig. [TH shows that, in the same 
way as for the experimental results, the conductivity is 
symmetric for V > and V < and it is practically inde- 
pendent of the tunneling current, although a significant 
population of excited states is present. Figure [TS] shows 
the excited state population as a function of the junc- 



tion bias for a conductivity of 2.0 10 ~ 6 S at V = and 
for an unpolarized tip. For large bias, the excited state 
populations are large (their sum is larger than the pop- 
ulation of the ground state); they are equal for positive 
and negative biases. Though no effect of these excited 
state populations is apparent in Fig. 1141 This is a direct 
consequence of Fig. [S] where the conductances of all the 
magnetic states are equal above the excitation thresholds 
in the case of unpolarised electrons (or holes) . Differences 
could only appear in the voltage region in between the 
inelastic thresholds, but in the present case, due to the 
smallness of the difference between the various excita- 
tion energy threshols, they are hidden by broadening ef- 
fects. All these qualitative features nicely agree with the 
experimental observations. The two sets though differ 
quantitatively, coherently with the decay rate compari- 
son in Fig. [6] The computed lifetimes are shorter and 
consequently, the effect of the excited state population 
on the junction conductivity appears for larger conduc- 
tances and larger tunneling currents. 

The effect of the populations of excited states on the 
polarized conductance is much visible in Figs. [12] and 
1131 which show the conductance as a function of bias 
for a fixed tip position. They also strongly emphasize 
a dcpcndancc of the conductance on the conductuctivity 
at V = 0. One can stress that if all the curves in e.g. 
Fig. [12] are plotted as a function of the junction current 
intensity, then they look very much alike, the popula- 
tion of the excited states depending dominantly on the 
current intensity . In such a plot, the role of the bias 
voltage and conductivity at V = only appear in the 
excitation threshold region. 
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FIG. 15: Population ol the six magnetic states in the 
Mn/CuN/Cu(100) system as a function of the tip bias for 
a non-magnetic tip. The magnetic field is along the x-axis 
and equal to 3T. The conductance at zero bias is equal to 
2.10 -6 S. In contrast to Fig. 1111 there is no bias asymetry. 



E. Symmetry of the tunneling process 

A very important feature in our strong coupling treat- 
ment of magnetic transitions is the symmetry of the tun- 
neling state, i.e. which is the value of St that is domi- 
nating tunneling through the adsorbate (see section IV Al 
above and a more detailed account in Ref. [HI)- As 
shown by our DFT study in Ref. [ll| , the St = 3 symme- 
try is dominating in the case of Mn/CuN/Cu(lll). Dif- 
ferent values of St lead to different strengths of the mag- 
netic transitions. For illustrative purpose, we checked 
the effect of changing the symmetry of the tunneling in- 
termediate in the present study. Figure [TO] shows the 
conductivity as a function of the STM bias, V, for a B 
field of 3T, a tip polarisation of i] = 0.24 and a conduc- 
tance at V = of 2.10 -6 S for three different tunneling 
symmetry hypothesis: i) St = 3 dominating, ii) St = 2 
dominating and Hi) both St = 3 and 2 contributing in 
a statistical way. This change of symmetry only con- 
cerns the transitions induced by the tunneling electrons 
and not the excited state decay by electron-hole pair cre- 
ation. Not surprisingly, the present physical situation 
with a partially polarised tip is highly sensitive to the 
symmetry of the tunneling intermediate. The three dif- 
ferent symmetries are seen to lead to different qualitative 
behaviours: the heights of the inelastic steps at small |V| 
are changing with the symmetry of the tunneling inter- 
mediate indicating a different efficiency of the various St 
symmetries in the magnetic excitation process. In addi- 
tion, depending on the symmetry and on the sign of V, 
the effect of excited state populations is seen to lead to 
an increase, a decrease or to a flat behaviour of the con- 
ductivity. Only the results for St = 3 do exhibit the 
qualitative behaviour found experimentally. This com- 
parison strongly supports our assignment of the St = 3 
symmetry as the dominant symmetry for the tunneling 
process in the present system. It also stresses the im- 
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FIG. 16: Relative conductance of the Mn/CuN/Cu(100) sys- 
tem as a function of the tip bias for various hypothesis in 
the spin-coupling scheme for the tunneling electrons (St — 3 
is the coupling scheme predicted by our DFT calculations). 
The tip polarization is r] = 0.24 and the B field, equal to 3 
T, is along the x-axis. The absolute conductance is equal to 
2.10 -6 S for zero bias. In the figure, the conductance is plot- 
ted in relative value, with the conductance for zero bias set 
to 1. The finite population of the excited states is taken into 
account. 



portance of considering the proper symmetry and proper 
spin coupling of the tunneling process for accounting for 
inelastic magnetic tunneling, especially with the present 
experimental protocol that involves a spin polarised tip 
and is thus much more selective. 



VI. CONCLUDING SUMMARY 

We reported on a theoretical study of the lifetime of 
excited magnetic states on single Mn adsorbates on a 
CuN/Cu(100) surface. Electrons tunneling through a 
single adsorbate are very efficient in inducing magnetic 
transitions as was revealed by IETS experiments; simi- 
larly, substrate electrons colliding on adsorbates are also 
very efficient in quenching a magnetic excitation in a sin- 
gle adsorbate on a surface, the magnetic excitation en- 
ergy being transferred into a substrate electron-hole pair. 
The present study on the quenching process is inspired 
by a recent study of magnetic excitations induced by tun- 
neling electron s 13 ' 14 . It is based on a DFT calculation of 
the system associated to a strong coupling approach of 
the spin transitions. We show that the decay rate of an 
excited magnetic state can be written as the product of 
the flux of substrate electrons hitting the adsorbate by 
a spin transition probability. In the Mn/CuN/Cu(100) 
system, the lifetime of the excited states is found to be 
rather short, the calculated lifetimes are typically in the 
0.04-0.4 ns range for an applied magnetic field in the 0-7 
T range. The lifetime decreases when the B field rises. 

For the comparison with Loth et al. experiments^, we 
performed a detailed study of electron tunneling through 
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the adsorbed Mn atom in the case of a polarized tip. 
Using polarized electrons (or holes) unveils several phe- 
nomena among which we can stress: 

• As shown in our strong coupling approach, the se- 
lection introduced by the symmetry of the total 
spin state in the tunneling process (total spin = 
spin of the tunneling electron + spin of the ad- 
sorbate) is very important and this leads to large 
differences in conductance between various situa- 
tions: tunneling electrons (or holes) with spin up 
(or down), adsorbates in different magnetic states. 
All these can be related to the existence of a domi- 
nant spin symmetry in the tunneling process, which 
is S T = 3 in the case of Mn on CuN/Cu(100) as 
shown here. 

• As the excitation probability by tunneling electrons 
is large, significant populations of excited states can 
exist for a finite current. Very interestingly, this 
leads to observable changes in the conductance as 
a function of current in the case of a polarized tip 
but not in the case of a non-polarized tip. This is 
a direct consequence of one of the key properties 
of the magnetic excitation process, as described in 
the strong coupling approach. The excitation pro- 
cess appears as a sharing process: a global tun- 
neling current is shared among the various mag- 
netic states, using sharing probabilities obtained 
from spin coupling coefficients. If we sum over all 
possible spin directions for the tunneling electron, 



these spin coupling coefficients sum to one. Thus, 
with a non-polarized tip, the conductance for biases 
above all inelastic thresholds is independent of the 
initial state and then independent of a finite sta- 
tionary population of excited states. In contrast, 
in the case of a polarized tip, the sum over electron 
spin directions is not complete and the conductance 
depends significantly on the magnetic state. 

• The present theoretical results reproduce all the be- 
haviours observed experimentally, basically those 
are consequences of the two phenomena described 
above. Quantitatively, although the excitation pro- 
cess is perfectly accounted for, the computed life- 
times of the excited states appear to be a factor 3 
shorter than those extracted from experiment. As 
a consequence, the experimentally observed vari- 
ation of the junction conductivity with the tun- 
nelling current is reproduced but for larger tun- 
nelling currents. 
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